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The notion of magnetic symmetry is reexamined in light of the recent observation of long range 
magnetic order in icosahedral quasicrystals [Charrier et al., Phys. Rev. Lett. 78, 4637 (1997)]. 
The relation between the symmetry of a magnetically-ordered (periodic or quasiperiodic) crystal, 
given in terms of a "spin space group," and its neutron diffraction diagram is established. In doing 
so, an outline of a symmetry classification scheme for magnetically ordered quasiperiodic crystals 
is provided. Predictions are given for the expected diffraction patterns of magnetically ordered 
icosahedral crystals, provided their symmetry is well described by icosahedral spin space groups. 
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Quasicrystals, which today are known to exist in 
perfectly-ordered thermodynamically-stable phases, have 
inspired a reexamination of the basic notions of crys- 
tallinity, long-range order, and symmetry [jl[. The recent 
observation by Charrier, Ouladdiaf, and Schmitt Q of 
long-range magnetic order in icosahedral quasicrystals of 
composition _RgMg42Zn5o (R — Tb, Dy, Ho, Er) is inspir- 
ing a similar examination of the nature of magnetic order 
which, in quasicrystals, had existed previously only as a 
theoretical construction || . The purpose of this Letter is 
to provide the first steps in this direction by explaining 
the notion of magnetic symmetry in a manner which ap- 
plies both to periodic and quasiperiodic crystals thereby 
establishing a tool for future analysis of neutron diffrac- 
tion data and providing a tentative interpretation for the 
findings of Charrier et al. ||. 

We choose to describe a magnetic material by its spin 
density field S(r). This field is a 3-component real- valued 
function, transforming like an axial vector under 0(3) 
and changing sign under time inversion. One may think 
of this function as defining a set of classical magnetic mo- 
ments, or spins, on the atomic sites of the material. The 
standard expression (as shown, for example, by Izyumov 
et al. p|) for the magnetic contribution to the intensity of 
elastic scattering of unpolarized neutrons is then given, 
in terms of the spin-density Fourier coefficients S(k), as 

/(k)(x|s(k)| 2 -|k.s(k)| 2 , (i) 

where k is the scattering wave vector and k is a unit 
vector in its direction. 

In making this choice we follow Litvin and Ope- 
chowski ^ who have developed a theory of "spin space 
groups" to describe the symmetry of periodic mag- 
netic crystals in terms of their spin density fields. In 
their theory, symmetry operations are those leaving 
the magnetic crystal invariant. These are the usual 
3-dimensional space-group operations — translations and 
proper or improper rotations — combined with rotations 
in "spin space" and possibly also time inversion. We ex- 
tend their theory here to deal with quasiperiodic crystals 
which possess neither translational symmetry nor, in gen- 
eral, any rotations that leave them invariant. We choose 



to extend the theory of spin space groups rather than two 
other commonly used theories for describing magnetic 
symmetry — those of color symmetry |EJ and the theory 
of representations of ordinary space groups — because 
of its direct predictions regarding the outcome of neutron 
scattering experiments. 

What should we expect to see in a neutron diffrac- 
tion diagram of a single magnetically-ordered crystal? 
To answer this question we must first clarify what we 
mean by "crystal." The International Union of Crystal- 
lography J8| defines a crystal to be "... any solid with 
an essentially discrete diffraction diagram." To be more 
concrete we consider spin density fields with well defined 
Fourier transforms, 

S(r) = J2 S(k)e* k - r , (2) 

keL 

in which the set L contains at most a countable infinity 
of plane waves. In a real experiment, due to the finite 
resolution of the apparatus, only a finite number of peaks 
whose intensities I(k) are above a certain threshold will 
be observed, resulting in a discrete diffraction diagram. 

What more can we say about the set of diffraction 
peaks beyond their being essentially discrete? We have 
shown elsewhere || that if S(r) describes a physically 
stable magnetic crystal, i.e., one which minimizes a suit- 
able Gibbs free energy, then the wave vectors k at which 
S(k) 7^ are closed under addition and subtraction, with 
the only exception of peaks that are required by symme- 
try to vanish. This implies in practice that once enough 
peaks have been observed, additional peaks will appear 
with increased experimental resolution only at integral 
linear combinations of peaks that already exist. 

This leads us to define the set L in Eq. (||) to be the 
set of all integral linear combinations of the wave vec- 
tors k determined by the diffraction diagram. We call 
this set the (reciprocal) magnetic lattice of the crystal. 
The rank D of L is the smallest number of wave vectors 
needed to generate it by taking integral linear combina- 
tions. As described above, we expect that in most (if not 
all) experimentally observed magnetic crystals this num- 
ber will be finite. Periodic 3-dimensional crystals have 
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rank D = 3; aperiodic 3-dimensional crystals have a rank 
D > 3; icosahedral quasicrystals, for example, have rank 
6. The first indication of the symmetry of the magnetic 
crystal is given by the set of (proper or improper) rota- 
tions which, when applied to the origin of Fourier space, 
merely permute the wave vectors of the magnetic lattice. 
This set forms a subgroup of 0(3) called the lattice point 
group Gl (also called the holohedry). 

The definition of "lattice" here is the same as in non- 
magnetic crystals, except that the magnetic lattice sup- 
ports the Fourier transform of a vector function S(r) 
whereas the lattice of a non-magnetic crystal supports 
the Fourier transform of a scalar function — that of the 
electronic or nuclear density p(r) of the crystal. We can 
therefore adopt the symmetry classification scheme, used 
in the non-magnetic case, also for megnatic lattices. This 
classification arranges lattices that have the same rank 
and the same lattice point group Gl into distinct Bra- 
vais classes according to the way in which their vectors 
transform under Gl (for details see, for example, Drager 
and Mermin JlO|]). Icosahedral lattices, for example, are 
arranged into three Bravais classes : P- lattices (prim- 
itive) contain all integral linear combinations of the six 
vectors 

= (±l,r,0); v( 2 < 5 > = (r,0,±l); V ( 3 < 6 > = (0,±l,r); 

(3) 

where r is the golden mean; F*-lattices (face centered in 
Fourier space) contain only those combinations in which 
the sum of the six integers is even; and /""-lattices (body 
centered in Fourier space) contain only those combina- 
tions in which the six integers are either all even or all 
odd. 

To say anything further about the nature of the 
diffraction peaks we must examine the symmetry of the 
spin density itself. We reformulate the theory of spin 
space groups by following the ideas of Rokhsar, Wright, 
and Mermin's "Fourier-space approach" to crystallogra- 
phy [fL2[ . At the heart of this approach is a redefini- 
tion of the concept of 3-dimcnsional point-group symme- 
try which enables one to treat quasicrystals directly in 
3-dimensional space The key to redefining point- 

group symmetry is the observation that certain rotations 
(proper or improper), when applied to a quasiperiodic 
crystal, even though they do not leave the crystal in- 
variant, take it into one that contains the same spatial 
distributions of bounded structures of arbitrary size. One 
finds that any bounded region in the unrotated crystal 
is reproduced some distance away in the rotated crystal, 
but there is, in general, no single translation that brings 
the two crystals into perfect coincidence. 

This weaker notion of symmetry, termed "indistin- 
guishability," is captured by requiring that any symme- 
try operation of the magnetic crystal leave invariant all 
spatially-averaged autocorrelation functions of its spin 
density field S(r) for any order and for any choice of 
components, 



= lim — / drS ai (n - r) ■ • • S an (r„ - r). (4) 

V^oo V Jy 

We have proven elsewhere (Appendix of E3)) that an 
equivalent statement for the indistinguishability of any 
two quasiperiodic multicomponcnt fields, in particular 
two spin density fields S(r) and S'(r), is that their Fourier 
coefficients are related by 

S'(k) = e 2mx ^S(k), (5) 

where x is a real- valued linear function (modulo integers) 
on L called a gauge function. Only in the case of periodic 
crystals can one replace 27rx(k) by k • d, reducing indis- 
tinguishability to the requirement that the two crystals 
differ at most by a translation d. 

We can now define the point group G of the magnetic 
crystal to be the set of operations g from 0(3) that leave 
it indistinguishable to within rotations 7 in spin space, 
possibly combined with time inversion. Accordingly, for 
every pair (5,7) there exists a gauge function, ^(k), 
called a phase function, which satisfies 

S( 5 k) = e 2 ™*3 (k) 7 S(k). (6) 

Since S([g/i]k) = S(g[/ik]), one easily establishes that 
the transformations 7 in spin space form a group T and 
that the pairs {g, 7) satisfying the point-group condition 
(^) form a subgroup of G x T which we call the spin 
point group Gs- The corresponding phase functions, one 
for each pair in Gs, must satisfy the group compatibility 
condition, 

VCff )7 ),(M)eG s : $^(k)^(/ik)+$^(k), (7) 

where "=" denotes equality modulo integers. A spin 
space group, describing the symmetry of a magnetic crys- 
tal, whether periodic or aperiodic, is thus given by a mag- 
netic lattice L, a spin point group Gs, and a set of phase 
functions $^(k), satisfying the group compatibility con- 
dition (0). 

In order to identify further the common symmetry 
properties of different magnetic structures, whose lat- 
tices and spin point groups are equivalent, one classifies 
their spin space groups into properly chosen equivalence 
classes called spin space-group types. This is achieved by 
organizing sets of phase functions satisfying the group 
compatibility condition (Q) into equivalence classes. Two 
such sets $ and 3>' are equivalent if: (I) they describe in- 
distinguishable spin density fields, related as in Eq. (|5|) 
by a gauge function %; or (II) they correspond to alterna- 
tive descriptions of the same crystal that differ by their 
choices of absolute length scales and spatial orientations. 
In case (I) $ and <!>' are related by a gauge transforma- 
tion: 

V( 5 , 7 )eGs: $'J(k)=^(k)+ X ( 5 k-k). (8) 
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TABLE I. Icosahedral lattice spin groups and their effect on the outcome of neutron scattering experiments. Nontrivial 
lattice spin groups are possible only with primitive (P) magnetic (reciprocal) lattices or with body centered (/*) magnetic 
(reciprocal) lattices. The lattice spin groups, which are groups of rotations in "spin space," are specified in the leftmost column 
in terms of their generating rotations, where 2^ is a 2-fold rotation about the i-axis in "spin space" (i = x,y, or z), 1' is the 
time inversion operation which takes every S(k) into — S(k), and 2'i is the product of the two. The form of S(k) issiven for 
each scattering wave vector in the magnetic lattice according to its indexing by the six vectors v' 8 ' defined in Eq. (p|). In the 
case of /"-lattices, three different scattering patterns are possible for each lattice spin group depending on the scale chosen 
for the generating vectors v' ! '. Entries in brackets are related to the ones above them through a scaling of the /'-lattice by 
the golden mean. The forms of the S(k) given in the table may be used in conjunction with expressions such as Eq. (|l|) to 
determine the outcome of neutron scattering experiments. 
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For a more rigorous definition of these equivalence crite- 
ria see the analogous classification of color space groups 
(Sec. Ill of Q). 

We said earlier that every wave vector k in the mag- 
netic lattice L of a magnetic crystal is a candidate for a 
diffraction peak unless symmetry forbids it. We are now 
in a position to understand how this happens. Given a 
wave vector k £ L we examine all spin point-group op- 
erations (g,"f) for which <?k = k. These elements form a 
subgroup of the spin point group which we call the lit- 
tle spin group o/k, Go. For elements (17,7) of Gg, the 
point-group condition (0) can be rewritten as 

7S(k) = e^WS(k). (9) 

This implies that every Fourier coefficient S(k) is re- 
quired to be a simultaneous eigenvector of all spin trans- 
formations 7 in the little spin group of k, with the eigen- 
values given by the corresponding phase functions. If a 
non-trivial 3-dimensional vector satisfying Eq. (^) does 
not exist then S(k) will necessarily vanish. It should be 
noted that the phase values in Eq. ([)]), are independent 
of the choice of gauge (||), and are therefore uniquely 
determined by the spin space-group type of the crystal. 

The process of determining the form of the simulta- 
neous eigenvector S(k) is greatly simplified if one makes 



the following observation. Due to the group compati- 
bility condition (^) the set of eigenvalues in Eq. (||) for 
all the elements (5,7) £ G| form a 1-dimensional repre- 
sentation of that group. Spin space-group symmetry thus 
requires the Fourier coefficient S(k) to transform under a 
particular 1-dimensional representation of the spin trans- 
formations in the little spin group of k. We also indepen- 
dently know that S(k) transforms under spin rotations 
as a 3-dimensional axial vector, changing its sign under 
time inversion. We therefore need to check whether the 
particular 1-dimensional representation, dictated by the 
spin space group, is contained within the 3-dimcnsional 
axial- vector representation. If it is not, then S(k) must 
vanish; if it is, then S(k) must lie in the subspace of spin 
space transforming under this 1-dimcnsional representa- 
tion. 

Of particular interest are spin transformations 7 that 
leave the spin density field indistinguishable without re- 
quiring any rotation in physical space. These transforma- 
tions are paired in the spin point group with the identity 
rotation e and form a subgroup of T called the lattice spin- 
group r e . The lattice spin group plays a key role in de- 
termining the outcome of elastic neutron scattering, for 
if a magnetic crystal has a nontrivial lattice spin group 
T e then {e} xf e C G| for every k in the magnetic lat- 
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tice, restricting the form of the corresponding S(k). In 
the special case of periodic crystals, the elements of T e 
are spin transformations that when combined with trans- 
lations leave the magnetic crystal invariant. The phase 
functions $2 (k) therefore contain the information which 
generalizes to the quasiperiodic case the so-called "spin 
translation groups" of Litvin and Opechowski || . 

What are the possible values of the phase functions 
$2(k)? Consider, for example, a lattice spin group gen- 
erated by an n-fold rotation 7 about the z-axis in spin 
space, with n > 2. Repeated applications of the group 
compatibility condition (0) to (e, 7)" = (e,e), where e is 
the identity in spin space, give 0=&2 s -)= n ^20 s -)- Thus 
$2(k)=i/n for some integer j. One can then easily verify 
through Eq. (|) that 

f(0,0,S 2 ) $2(k)=0, 
S(k)= (S x ,±iS x ,0) $2(k)=±i (10) 
I (0, 0, 0) otherwise. 

We have enumerated the distinct lattice spin groups 
for icosahedral quasicrystals and have found that non- 
trivial lattice spin groups are possible only with P- or 
/■"-magnetic lattices. For each of the lattice spin groups 
we have also determined the expected form of S(k), that 
is required through (||) by the spin space group sym- 
metry, for every scattering wave vector k in the mag- 
netic lattice. These results are summarized in Table |. 
Further restrictions (not tabulated here) may exist for 
wave vectors k lying in the invariant subspaces of non- 
trivial point-group operations. We emphasize that the 
diffraction patterns are described here as magnetic lat- 
tices with missing points rather than nuclear lattices that 
are shifted by so-called "magnetic propagation vectors." 

Charrier et al. [Q observe magnetic reflections for the 
i?sMg42Zn5o quasicrystals at wave vectors of the form 
k = y^ -_ i mjV W + pv^\ where the vW are defined in 
Eq. (y), the mi are either all even or all odd, j — 1, . . . , 6, 
and p — ±i. The nuclear reflections form a body- 
centered icosahedral (reciprocal) lattice, obtained from 
the expression above but with p = 0, corresponding to 
face-centered ordering in direct space. If the magnetic 
structure indeed has icosahedral symmetry, and is not 
merely a collection of magnetic domains in which icosa- 
hedral symmetry is broken, then the magnetic lattice, 
which is formed by taking all integral linear combina- 
tions of the observed magnetic reflections, is a primi- 
tive icosahedral lattice containing all vectors of the form 
k = Yli=i n i(f v ^) with an Y integers rij. All reflections 
at wave vectors with Y^=i n i even are not observed by 
Charrier et al. which is consistent with having a lattice 
spin group 1', as shown by the first entry for P-lattices 
in Table Q. Other peaks that are not observed might just 
be too weak in the current experiment rather than ac- 
tually missing due to symmetry requirements. This may 
suggest that in direct space the magnetic structure has 
an underlying antiferromagnetic body-centered icosahe- 
dral ordering (analogous to the chemical ordering in the 



cubic cesium chloride structure). This may occur, for 
example, if only a fraction of the rare-earth atoms have 
their magnetic moments aligned, this fraction arranged 
in a body-centered icosahedral super-structure made of 
tiles with twice the edge length. One would need to de- 
velop actual models to test such hypotheses. 

We have demonstrated that through its fairly simple 
selection rules (||) the theory of spin space groups pro- 
vides a valuable tool for analyzing neutron diffraction di- 
agrams of either periodic or quasiperiodic magnetically- 
ordered crystals. This is not to say that the use of color 
symmetry, which has been extended to quasicrystals [jlJJ , 
and the use of representations of ordinary space groups, 
which have yet to be extended to quasicrystals, will not 
offer any additional insight. 

This paper is dedicated to Hans-Ude Nissen on the 
occasion of his 65th birthday. The author thanks Michael 
Cross, Veit Elser, Peter Weichman, Benoit Charrier, and 
Denys Schmitt for valuable discussions, and the latter 
two also for sending their preprints prior to publication. 
This work was supported by Caltech through a Division 
Research Fellowship in Theoretical Physics. 
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